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Abstract. The LINER nucleus of the nearby spiral 
galaxy M81 was pointed by BeppoSAX, which caught it 
at the highest (2-10) keV flux level observed so far. The 
LEGS, MEGS and PDS data, extending over (0.1-100) 
keV, are used to investigate the physical similarities and 
differences between LINERs and AGNs. The continuum 
is well fitted by a power law of photon index T ^ 1.84, 
modified by little absorption due to cold material; this 
extends to (0.1-100) keV the validity of a similar ASCA 
result. Superimposed on the continuum BeppoSAX de- 
tects a 6.7 keV emission line (confirming another ASCA 
result) and an absorption edge at ^ 8.6 keV. Both spec- 
tral features are consistent with being produced by iron 
at the same high ionization level, and probably also with 
the same column density. So, we suggest that they origi- 
nate from transmission through highly ionized thin mate- 
rial. Goncerning the origin of the continuum emission, we 
do not observe signs of reflection from the optically thick 
material of an accretion disk, as usually found in Seyfert 
I's (a 6.4 keV emission line and a broad bump peaking 
at 10-20 keV). The low bolometric luminosity of the nu- 
cleus of M81 is consistent with being produced by advec- 
tion dominated accretion; in this case the X-ray emission 
should be dominated by Gomptonization, rather than by 
bremsstrahlung, in order to reproduce the steep spectrum 
observed over the (0.1-100) keV band. 
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1. Introduction 

M81 (NGG3031) is a nearby spiral galaxy with a promi- 
nent bulge and well defined spiral arms, morphologically 
similar to M31 (Table 1). It has been well studied at all 
frequencies, from radio (e.g., Beck et al. 1985, Bietenholz 
et al. 1996), to optical and UV (e.g., Ho et al. 1996), to 
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X-rays (Sect. 2). On the basis of these observations, M81 
turned out to be the closest galaxy to show the spectro- 
scopic signatures of a LINER (Ho et al. 1997). In addi- 
tion, some observed properties make the nucleus of M81 
a good candidate for a low luminosity AGN (LLAGN). 
These include the presence of a broad component of the 
Ha emission line (Peimbert & Torres-Peimbert 1981), a 
compact radio core (Bietenholz et al. 1996), a pointlike X- 
ray source coincident with the optical nucleus (Fabbiano 
1988, hereafter F88), and the presence of a power law con- 
tinuum in the 2-10 keV energy band (Ishisaki et al. 1996, 
hereafter 196). Dynamical studies also suggest the pres- 
ence of a super-massive object at the galaxy nucleus, of 
mass 4 x IO^Mq (Ho 1999). 

Here we report the results of an analysis of the prop- 
erties of the nucleus of M81 over 0.1-100 keV, and discuss 
them in relation to previous X-ray observations (mostly in 
the 2-10 keV band) and the unsolved issue of sorting out 
the physical similarities and differences between LINERs 
and AGNs. In particular the origin of LINERs is still un- 
certain: their emission line spectrum could be powered by 
a LLAGN, in which case they could represent the miss- 
ing link between normal galaxies and the less luminous 
AGNs, or by starburst activity (Terlevich et al. 1992). The 
very large energy band of BeppoSAX coupled to its spec- 
tral resolution (8% FWHM at 6 keV) is especially suited 
for reliably measuring the high energy continuum, and to 
search for a thermal component and Fe-K features. There- 
fore we can investigate more accurately than with previous 
instruments the similarities and differencies of the X-ray 
properties of LINERs with those of AGNs. 



2. X-ray observations of M81 prior to BeppoSAX 

M81 has been pointed by many X-ray satellites. It was first 
observed in X-rays with Einstein (Elvis & Van Speybroeck 
1982; F88). This imaging satellite revealed the presence of 
several discrete sources in the M81 region, the brightest 
of which coincides in position with the nucleus (Fig. 1). 
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Table 1. General characteristics of M81 



Type"^ 


RA 


Dec 




dO a 


logl/B 


iVH,Gal " 




(J2000) 


(J2000) 


(Mpc) 


(mag) 




(cm-2) 


Sab 


9''55™33':2 


69°03'55" 


3.6 


7.39 


10.32 


4.1x10^° 



" from de Vaucouleurs et al. (1991). is the total B magnitude, corrected for Galactic and internal extinction. 
Cepheid-based distance (Freedman et al. 1994). 

total B-band luminosity, derived using the indicated distance and B^. 
Galactic neutral hydrogen column density from Stark et al. (1992). 



Its spectral representation was given in terms of a steep 
power law (F = 4.0) or in terms of a thermal component 
of kT = 1.1 keV, both absorbed by a hydrogen column 
density A'h in excess of the Galactic value A^H.Gai (in Ta- 
ble 1). Subsequent observations by GINGA (Ohashi et al. 
1992), BBXRT (Petre et al. 1993) and ROSAT (Boiler 
et al. 1992) gave a spectral interpretation in terms of a 
power law with F = 2 and Nu ranging from being con- 
sistent with A^H.Gai to ten times higher. Petre et al. also 
found a thermal component of kT — 0.4 keV. 

ASCA pointed M81 several times between 1993 and 
1999 (196, lyomoto 1999). It was found that, over such a 
long period, the 2-10 keV luminosity of the nucleus was 
highly variable (Fig. 2). However no spectral variability 
was found, and the average spectrum was well fitted by a 
power law continuum of photon index F ^ 1.85, absorbed 
by a column density of iVn ^ 10^^ cm^^, plus a thermal 
component with a temperature of kT ^ 0.6 — 0.8 keV. A 
possibly broad or complex iron emission line centered at 
6.6-6.9 keV, with an equivalent width of 170!!igQ eV, was 
also detected. 



3. X-ray data analysis 

M81 was observed by three of the Narrow Field In- 
struments onboard the BeppoSAX satellite (Boella et 
al. 1997a): the Low Energy Goncentrator Spectrometer 
(LEGS), the Medium Energy Goncentrator Spectrometer 
(MEGS), and the Phoswich Detector system (PDS; Fron- 
tera et al. 1997). The journal of the observation is given 
in Table 2. The LEGS and MEGS are grazing incidence 
telescopes with position sensitive gas scintillation propor- 
tional counters in their focal planes. The MEGS, which 
at the epoch of the observation consisted of two equal 
units, has a field of view of 56' diameter, and works in the 
range 1.6-10 keV (Boella et al. 1997b). The LEGS oper- 
ates at softer energies (0.1-4.5 keV), has a field of view of 
37' diameter, an energy resolution a factor of ^ 2.4 bet- 
ter than that of the ROSAT PSPC, but an effective area 
much lower (between a factor of 6 and 2 lower, going from 
0.3 to 1.5 keV; Parmar et al. 1997). The PDS is a coUi- 
mated instrument, operating in rocking mode (i.e., half of 
the time on the source and half on the background direc- 



tion), that covers the 13-300 keV energy band. It has a 
triangular response with FWHM of ~ 1°3 (Frontera et al. 
1997). 

The cleaned and linearized data have been retrieved 
from the BeppoSAX Science Data Genter archive, and 
later reduced and analysed using the standard software 
(XSELEGT vl.4, FTOOLS v4.2, IRAF-PROS v2.5, and 
XSPEG vlO.O). For the MEGS, we used the event file made 
by merging the data of the two properly equalized MEGS 
units. The PDS data reduction was performed using in- 
dependently both XAS (v. 2.0, Ghiappetti & Dal Fiume 
1997) and SAXDAS (v.1.3.0, Fiore et al. 1999) software 
packages, and yielded consistent results. In the following, 
we will refer to the results obtained with the SAXDAS 
package. 

3.1. Spatial analysis and extraction regions 

The satellite pointed the optical center of the galaxy (in 
Table 1). The center of the X-ray emission was found to 
be located at 9''55™26':5, 69°03'43" for the LEGS, and 
9''55"32^5, 69°03'35" for the MEGS, i.e., within - 1' from 
the optical center of the galaxy given in Table 1. There- 
fore, within the accuracy with which positions are given 
by BeppoSAX , the X-ray centers coincide with the M81 
nucleus. The LEGS and MEGS images in the 0.1-4.0 and 
1.65-10.5 keV energy bands respectively are shown in Fig. 
1. Also shown are the positions of the sources detected in 
the field by the Einstein observation. 

For a point source, the PSF of the MEGS in- 
cludes 80% of photons of energies > 1.5 keV within 
a radius of 2^7 (BoeUa et al. 1997b). The PSF of 
the LEGS is broader than that of the MEGS below 
1 keV, while it is similar to it above 2 keV (see 
bttp://www.sdc.asi.it/software/cookbook). In order to 



study the spectral properties of the nucleus of M81, given 
these instrumental characteristics, counts have been ex- 
tracted from a circle of 3' radius for the MEGS, and of 
4' radius for the LEGS. The larger radius adopted for the 
LEGS is motivated by its larger PSF at softer energies. 
A radius of 3' was also adopted for the extraction of the 
Einstein and ASCA spectra (F88, 196). We estimate in 
Sect. 3.6 the contribution of the unresolved emission from 
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Table 2. Observation Log 



Date 


Exposure time" (ks) 


Count Rate'' (ct/s) 






LEGS MEGS PDS 


LEGS MEGS 


PDS 


1998 Jun 4 


43.20 100.3 36.40 


0.329±0.003 0.496±0.002 


0.57±0.04 



° On-sourcc net exposure time. The LEGS exposure time is considerably shorter than the MEGS one, because the LEGS can 
operate only when the spacecraft is not illuminated by the Sun. 
Background subtracted source count rates, with photon counting statistics errors. 



Fig. 1. The LEGS and MEGS images of M81, smoothed with a gaussian of cr = 32". Gontours plotted are logarith- 
mically spaced from 0.7 to 90% of the peak intensity, for the LEGS, and from 0.6 to 90% for the MEGS. The labels 
X1-X12 mark the positions of the sources detected by Einstein in the field (F88); X5 is the nucleus of M81. The 
optical image of M81 is from the Digital Sky Survey of the Space Telescope Science Institute. 



the galaxy to the LEGS and MEGS spectra extracted as 
described above. 

The starburst galaxy M82, that lies - 37' from M81, 
is inside the PDS field of view (and outside the LEGS 
and MEGS ones). Given the PDS triangular response (1°.3 
FWHM), half of the M82 flux between 13-300 keV could 
"contaminate" the M81 observation. A BeppoSAX obser- 
vation of M82 performed one year earlier indicates, how- 
ever, that M82 has a steeper spectrum than M81 (photon 
index of 3.8 ±1.8 in the PDS energy band), and a 13-50 
keV flux <1.7 X 10"" erg cm^^ (Gappi et al. 1999), 
i.e., a factor of >2 lower than that of M81 in the same 
energy range. We estimate then that the PDS spectrum 
of M81 could be contaminated by at most 25% of its flux. 
Therefore, we added a systematic uncertainty of this level 
to the spectral data between 13-50 keV. 



3.2. Variahility 

The long term variation of the M81 (2 10) keV flux, ob- 
tained by collecting observations from various observato- 
ries, is plotted in Fig. 2. During the BeppoSAX pointing 
the nucleus was at the highest level ever observed. From 
1993 to 1998 it brightened by a factor of ~ 4. 

A short term variation is clearly detected by 
BeppoSAX] the light curves obtained from the LEGS and 
MEGS data are plotted in Fig. 3. The count rates of both 
instruments show a variation of ^ 30% from peak to val- 
ley over roughly one day. The fit of the MEGS light curve 
with the sine function (Fig. 3, bottom) is significantly bet- 
ter than that with a constant value; the period turns out 
to be almost 2 days (43 hours). Since this period is com- 
parable to the length of the observing time, the presence 
of a periodical change cannot be established; future ob- 
servations might reveal whether this periodicity is real. 
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Fig. 2. Long term variation of the (2 10) keV flux from 
the nucleus of M81. The points refer to observations by 
Einstein, EXOSAT, GINGA, ROSAT, BBXRT, ASCA, 
and BeppoSAX. Typical 90% confidence errors on fluxes 
are not larger than ±10%. The fluxes coming from ASCA 
data have been kindly provided by N. lyomoto (1999), and 
that from EXOSAT data by P. Barr (priv. comm.). For 
the (2-10) keV fluxes of Einstein and ROSAT we have 
adopted the extrapolations by Petre et al. (1993). 

The BeppoSAX light curve extends over a time interval 
as long as the longest of the ASCA observations (made 
during Aprfl 1993; Serlemitsos et al. 1996). The ASCA 
light curve in that case showed variability of the order of 
~ 20% on a timescale of ~ 1 day. 

It is unlikely that the flux variation is produced by an 
X-ray binary in the field, because it should be extremely 
luminous (Lx ~ 2 x 10"'° erg s~^ in the 0.1-2 keV band, 
and Lx ~ 10*° erg over 2-10 keV). No sources so 
luminous have been detected by Einstein and ROSAT 
(see Sect. 3.6). Another fact points against the hypothesis 
of a binary responsible for the variability: the variation 
detected by ASCA, which had observed the nucleus of 
M81 during a phase of lower flux level in 1993 (Fig. 2), had 
a smaller amplitude than that detected by BeppoSAX; a 
variation of larger amplitude is instead expected when the 
flux of the nucleus is lower, if it is produced by a binary. 

3.3. Spectral analysis 

The background spectrum was estimated from blank fields 

event files (released in Nov. 1998). These blank-fields 
events were accumulated on five different pointings of 
empty fields; the extraction regions used for them corre- 
spond in size and position to those of the source. Spectral 
channels corresponding to energies 0.12-4 keV, 1.65-10.5 



MECS 




d 



5x10* 10^ 1.5x10^ 

Time (s) 



Fig. 3. Top: LEGS light curve; counts are from 4' radius, 
and 0.1-4 keV. Bottom: MECS light curve; counts are from 
3' radius, and 2 10 keV. The sine function that best fits 
the data is also plotted. 



keV, and 15-200 keV have been used for the analysis of 
the LEGS, MECS and PDS data respectively. 

The data have been compared to models, convolved 
with the instrumental and mirror responses, using the 
minimization method. For this comparison the original 
channels have been rebinned in order to sample the in- 
strument resolution with the same number of channels at 
all energies, and in order to have bins adequately filled for 
applicability of the statistic to assess the goodness of fit 
(Gash 1979). The spectral response matrices and effective 
area files released in September 1997 have been used in the 
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fitting process. Wc fitted the models simultaneously to the 
LEGS, MEGS and PDS data. In the fitting two normal- 
ization constants have been introduced to allow for known 
difi'erences in the absolute cross-calibrations between the 
detectors (Fiore et al. 1999). The models are modified by 
photoelectric absorption of the X-ray photons due to inter- 
vening cold gas along the line of sight, of column density 
VVh- The results of the spectral analysis are presented in 
Table 3. 

3.3.1. Fitting results 

A single thermal model cannot reproduce the broad band 
data, in which the nucleus is detected up to 100 keV. The 
LEGS and MEGS data constrain the temperature to be 
around 8 keV, but with this temperature the model of 
course totally fails to reproduce the PDS data. On the 
other hand, a thermal component that fits the high energy 
region would have kT ^ 90 keV, and would be too flat to 
reproduce the LEGS and MEGS data. A two-component 
thermal model (with temperatures kT ~ 3 and kT ^ 60 
keV) gives a better fit, but still of poor quality. A single 
power law component of photon index F ~ 1.85 can in- 
stead reproduce the data over the whole energy band (Fig. 

4) . As can be seen looking at the residuals (see also Fig. 

5) , excess emission is present between 6-7 keV, as already 
found by ASCA (Sect. 2). Adding a narrow gaussian line 
to the model gives a decrease of of 21 for two additional 
parameters (the line center energy and normalization); so, 
this line is statistically significant at > 99.9% confidence 
level in an F-test (this follows the treatment, specific for 
fitting of X-ray spectral data, about how to assess the use- 
fulness of adding extra parameters to an initial model by 
Mahna et al. 1976). The fine is centered at 6.70to keV, 
and its equivalent width is EW=104l4Q eV (see model A 
in Table 3; Fig. 7). If the width a of the line is left as 
a free parameter, wc obtain for it at the best fit a value 
of 110 eV and an upper limit at 90% confidence for one 
interesting parameter of 300 eV. The 6.4 keV value for the 
line center energy is excluded at > 99% confidence level. 
Its introduction in the fitting provides a 90% confidence 
upper limit on its EW of 42 eV, and no improvement of 
the fit. 

The power law is absorbed by cold material exceed- 
ing the Galactic column density by Nji ~ 7 x 10^° cm~^. 
This intrinsic column density is consistent with that de- 
rived optically: Filippenko & Sargent (1988) estimated an 
E(B-V)=0.ltg i^ for the broad line region of M81, which 
translates into A^n = 5.8 x 10'^'^ E{B - V) = 5.8 x lO^o 
cm~^. 

A further improvement in the quality of the fit is ob- 
tained by fitting with a power law plus a thermal com- 
ponent (model B in Table 3; the mekal model describes 
the thermal emission from an optically thin hot plasma, 

both from contimmm and lines). Then F ^ 1.85 and again 
Nn ~ 10^^ cm^^ for the power law component; for the 



Table 3. Spectral analysis 



Model parameters 




Best fit values 


A. pow-|-gauss: 






10-2° iV5 (cm-2) 




7.2 (6.1-8.5) 


r 




1.86 (1.84-1.89) 


10"F(crg cm-2 s) 




2.1, 3.8, 7.4 


Line E(keV) 




6.70 (6.59-6.81) 


Line EW (eV) 




104 (64-133) 


A. 1 ^ 




193/178 


B. mekal^-hpow: 






Z (Zo) 




0.5 (0.4-1.3) 


kT (keV) 




6.5 (3.5-7.9) 


10"F^efcai(erg cm" 


s) 


0.4, 0.7, 0.6 


10-2° ATg (cm-2) 




7.1 (5.7-8.2) 


r 




1.85 (1.76-1.88) 


10"Fpc.„,(erg cm-2 




1.7, 3.1, 5.6 






181/177 


C. pow-|-gauss-|-edg 


:e: 




10-2° ATg (cm-2) 




6.7 (5.7-7.9) 


F 




1.84 (1.82-1.87) 


10"F(crg cm-2 s) 




2.1, 3.8, 7.4 


Lino E(kcV) 




6.69 (6.58 6.83) 


Line EW (cV) 




84 (48-121) 


Edge E(koV) 




8.6 (7.8-9.0) 






0.15 (0.07-0.24) 


xV^ 




182/176 


D. mekal''-|-pow-|-gauss-|-edge: 








0.03 (0.-0.09) 


kT (keV) 




0.50 (0.25-0.95) 


10^^Fmefcai(erg cm' 


' s) 


0.3, 0, 


10-20 ATg (cm-2) 




12.0 (9.0-15.6) 


r 




1.86 (1.83-1.89) 


10"Fpo„(erg cm-2 


s) 


1.8, 3.8, 7.2 


Lino E(kcV) 




6.70 (6.58-6.83) 


Line EW (eV) 




95 (49-123) 


Edge E(keV) 




8.6 (7.8-9.1) 






0.14 (0.06-0.23) 






172/173 



" Column density of neutral hydrogen in addition to A^n.Gai- 

The column density of neutral hydrogen absorbing the mekal 
component is fixed at the Galactic value of 4.1 x W^''cm~^. 

Note: v is the number of degrees of freedom of the fit. The 
values between parenthesis, close to the best fit values, give the 
90% confidence intervals for one interesting parameter (these 
intervals correspond to Ax^ ~ 2.71). 

Note: three values for the absorbed fluxes F are given, respec- 
tively for the (0.1-2), (2-10), and (10-100) keV energy bands. 

thermal component kT ^ 6 keV, A''h ~ iVn.Gai (imposed, 
otherwise it would be lower) and Z ~ 0.5Zq. This model 
gives a fit of better quality than model A because the 
thermal component reproduces the 6.7 keV line and gives 
slightly smaller residuals than a simple power law in a 
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power law 




1 10 100 

channel energy (keV) 

Fig. 4. BeppoSAX LEGS, MEGS and PDS observed spectra of M81 (crosses), modeled with a power law of T = 1.85 
(solid line). The residuals, in terms of a's, between the data and the model are plotted below. 



channel energy (keV) 



Fig. 5. The data/model ratio for the best fit continuum 

model (the power law shown in Fig. 4) around the 6.7 keV 
emission line and the 8.6 keV absorption edge. 



few channels below 1 kcV. This thermal component con- 
tributes 16% of the total unabsorbed flux over 0.1-10 keV. 



3.3.2. Detection of an ionized absorber 

There is another spectral region, located around 9 keV, 
that shows residuals when fitted with a simple power law 
(Figs. 4 and 5). In fact the quality of the "power law plus 
gaussian line" fit improves by adding an absorption edge 
(Fig. 6; model C in Table 3). This turns out to be located 
&tE = S.eto.t keV, with optical depth r = 0.151^;^^ (Fig. 
8) . The presence of one edge is statistically significant at 
> 99% confidence level (from an F-tcst; decreases by 11 
for 2 more fitting parameters). We cannot exclude that the 
'hole' around 9 keV could be produced by more than one 
ionization stage. Fitting with two edges, at energies fixed 
at 8.5 keV (FeXXIII) and 8.8 keV (FeXXV), improves 
the fit slightly but not significantly. This is the first time 
an absorption edge is detected in the X-ray spectra of 
LINERs. 



3.4- Spectral analysis and variability 

We have also looked for spectral differences between the 
two states of low and high flux, i.e., we have studied sep- 
arately the spectra obtained from the first and the last 
~ 7.5 X 10^ s of the pointing (see Fig. 3). We do not find 
significant changes between the two states in any spec- 
tral property described above (power law slope, emission 
line and absorption edge center energies, optical depth of 
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power law + gaussian line + absorption edge 




1 10 100 

channel energy (keV) 

Fig. 6. BeppoSAX LEGS, MEGS and PDS observed spectra of M81 (crosses), modeled with a power law plus gaussian 
line plus absorption edge (model G in Table 3). The residuals, in terms of ct's, between the data and the model are 
plotted below. 




6,5 6.6 6.7 6.B 6.9 7 7 7.5 B 8.5 9 9.5 

LineE keV sdgeE keV 



Fig. 7. The 68%, 90% and 99% confidence contours for Fig. 8. The 68%, 90% and 99% confidence intervals for 
two interesting parameters for the line energy and normal- two interesting parameters for the energy and optical 
ization (in units of photons cm~^ s~^ in the line). depth r of the absorption edge. 



the edge). Going from the low to the high state the nor- 
malization of the power law increases, and the equivalent 

width of the 6.7 keV line marginally decreases. More pre- 
cisely, when a fit with a power law -|- gaussian line -|- 



edge is performed, going from the low to the high state 
r slightly decreases from 1.86l° to 1.84^°;°^, while the 
EW drops from ISOtlt cV to bbt^ cV. The line cen- 
ter energy goes from 6.73 to 6.65 keV, with the associated 
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confidence intervals almost overlapping. The line flux goes 

from 4.4 x 10~^ pliotons cm^^ s^^ to 2.1 x 10~^ photons 

— ? — 1 
cm s . 

3.5. X-ray fluxes and luminosities 

In Table 3 we report the observed fluxes in the energy 
bands (0.1-2), (2-10) and (10-100) keV. The luminosities 
are obtained from fluxes corrected for the observed total 
A^'h, adopting a distance of 3.6 Mpc. So, for the model 
power law + gaussian line + edge (model C), in the three 
bands Lx = 7.5, 4.3, 9.4 x 10^" erg s'^. The (0.1-100) 
keV luminosity is 2.1 x 10*^ erg . 

3. 6. Background emission from the M81 galaxy 

Before ascribing the derived spectral properties to the nu- 
cleus of M81, we must evaluate the contribution, within 
the region used to extract the LEGS and MEGS spectra, 
from the various possible non-nuclear X-ray components 
of M81: stars, supernova remnants, X-ray binaries, and 
possibly a hot ISM. Einstein resolved a few sources to 
within 3' (Fig. 1). In particular, we expect a contribution 
from sources X2, X3, X6, and X7. Given their positions, 
their Einstein fluxes (F88), and the shape of the PSF, we 
estimate their contribution to be < 6% of the total net 
counts within a radius of 3' from the nucleus^. The point 
sources resolved in the ROSAT HRI image within 3' of 
the nucleus are 11 (G. Fabbiano, private communication). 
Their HRI count rates give a contribution of <6% to the 
total (0.1-2) keV flux within a circle of 3' radius. Note 
that there are no sources bright enough to produce the 
short term variability shown by BeppoSAX (Fig. 3, Sect. 
3.2). From an estimate kindly provided by G. Fabbiano of 
the 'galaxy background' contribution (which includes also 
the ROSAT unresolved emission) we derive that >85% 
of the total counts within 3' belong to the active nucleus, 
in the ROSAT (0.1-2) keV band; this calculation takes 
into account also the increase in flux observed between 
the ROSAT and BeppoSAX pointings (Fig. 2). 

The galaxy background contribution is difficult to esti- 
mate with data from satellites sensitive to higher energies 
than ROSAT, due to their poor angular resolution. This 
contribution has proven to be not influential for the spec- 
tral analysis in the ASCA (0.5-10) keV band (196). 196 
estimated that the galaxy background outside the central 
3' is well fitted by a bremsstrahlung of kT ^ 8 keV (which 
led 196 to ascribe the galaxy background to the integrated 
emission from LMXBs). A spectral fit of the nuclear re- 
gion with such a bremsstrahlung component restricted the 

^ During the Einstein pointing the source X5 was in a low 
state, by at least a factor of ten lower than what was found later 
by EXOSAT, GINGA, ROSAT, BBXRT, and a factor of ~ 
fourty lower than found by BeppoSAX (Fig. 2); so, using the 
Einstein fluxes for the various point sources we can estimate 
an upper limit to their contribution to the total net counts. 



galaxy background contribution to < 12% of the 2-10 keV 
luminosity of the active nucleus. A similar analysis for the 
LEGS and MEGS data would not give results more reli- 
able than those obtained with ASCA, due to the angular 
resolution of these instruments and the vignetting effect 
that becomes large off-axis. All in all, we are confident that 
> 80% of the (0.1-100) keV net counts extracted from a 
circle of 3' radius from the galaxy center come from the 
active nucleus. 

4. Discussion 

The simplest model that explains all the features of the 
0.1-100 keV spectrum is a power law of F '--^ 1.84, to which 
a gaussian line at 6.7 keV and an absorption edge at ~ 8.6 
keV are added. The energy of the line indicates Ka emis- 
sion from He-like iron (FeXXV), while its 90% confidence 
interval goes from FeXXII to FeXXV. The absorption edge 
can be produced by ions from FeXVII to FeXXV, within 
the 90% confidence interval, and its best fit energy cor- 
responds to the FeXXIV K-edge (^ 8.6 keV; Makishima 
1985). Both features therefore come from highly ionized 
material. Moreover they could be consistent with each 
other also in column density. The best fit value for the 
optical depth r, when assuming a cosmic abundance of 
iron, and the photoionization cross section for the K-shell 
for ionization stages from FeXXIII to FeXXV (Krolik & 
Kallman 1987), corresponds to a hydrogen column density 
of ~ 2 X 10^^ cm~^ for the ionized material. The observed 
EW of ^ 100 cV for the Fe-K line can be produced by 
reflection from or transmission through a column density 
of A^H ~ 10^'^ cm^^, in the case of neutral material (Mak- 
ishima 1985). In case of highly ionized material, an EW 
of ^ 100 eV is expected to correspond to slightly lower 
column densities than for the neutral case (Turner et al. 
1992). 

As can be seen from Table 3, a model made by the 
superposition of a power law and a thermal component 
at kT ~ 6 keV (model B) is equally statistically accept- 
able as a model with a power law -I- gaussian line -I- edge 
(model G), because it reproduces the 6.7 keV line, and 
the residuals in a few channels below 1 keV are slightly 
smaller. Since the absorption edge remains totally unex- 
plained with a thermal model, we prefer to adopt model 
G. To reduce its residuals below 1 keV, the addition of 
a soft {kT ^ 0.5 keV) thermal component proved effec- 
tive (model D in Table 3); it contributes ~ 12% of the 
(0.1-2) keV absorbed flux. This soft component is not re- 
quired by the data though, since is reduced by ten, 
but the number of free parameters increases by three. A 
soft component had been detected also in the analysis of 
BBXRT data by Petre et al. (1993), and by 196 (Sect. 
2). This soft thermal emission could represent the 'galaxy 
background' within 3' of the nucleus, that we estimate to 
account for <15% of the total emission within that ra- 
dius in the ROSAT band (Sect. 3.6). Soft emission with 
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this temperature and sub-solar abundance is common to 
spiral galaxies, and may well come from heterogeneous 
origins (supernova remnants, stars, thermal gas created 
in a starburst; e.g., Serlemitsos et al. 1996, lyomoto et al. 
1997). 

The possible origin of the adopted model (power law 
+ gaussian line + edge) is discussed in the following. The 
power law component that well reproduces the (0.1-100) 
keV spectrum is usually related to the presence of an 
AGN. This presence is not unexpected here: based on the 
velocities and dimensions of the broad line region, and as- 
suming that gravity dominates the motion of the clouds. 
Ho et al. (1996) derived a central mass of (0.7-3) x IO^Mq 
for the nucleus of M81, recently updated to 4 x IO^Mq 
(Ho 1999). Moreover, no absorption much in excess of the 
Galactic value is seen in the X-ray spectrum. This prop- 
erty, together with the value of F close to that found for 
Seyfert I's (Nandra & Pounds 1994), makes M81 similar 
to Seyfert I's. 



law continuum, at 0.74 and 0.87 keV; we do not observe 
these edges, because the material is highly ionized. Con- 
cerning the location of such material, it could be close to 
the source of ionizing photons. In the hypothesis that its 
line of sight thickness is much less than d, its distance from 
the source of photoionizing photons, Reynolds & Fabian 
(1995) derive a constraint on d: d < L/N^, where L is the 
luminosity of ionizing photons, N is the column density 
of photoionized matter, and ^ is its level of ionization^ 
(^ = L/nd'^, where n is the number density of photoion- 
ized matter, Kallman & McCray 1982). So d < lO^'^ cm, if 
~ 2 X 10^'^ cm~^ (as derived from the absorption edge 
in Sect. 4) and L ^ 2x 10^^ erg s~^ for the observed (0.1- 
100) keV luminosity of the nucleus of M81. This places 
the absorber slightly closer to the nucleus with respect to 
what is found typically for Seyfert I's, where the so-called 
warm absorber is located at radii coincident with, or just 
outside, the broad line region (e.g., Reynolds & Fabian 
1995). 



Origin oj the iron-K emission and absorption 

Are the 6.7 keV emission line and the 8.6 keV absorption 
edge produced by reflection (e.g., from an X-ray illumi- 
nated disk, Zycki & Czerny 1994), or by transmission? 
And where is the highly ionized material responsible for 
the line and the edge located? In the case of reflection 
from a disk, its ionization state should be higher than in 
Seyfert I's, because the iron line energy is higher (6.7 keV 
instead of 6.4 keV). But there are both observational and 
physical problems with this hypothesis. The observational 
problem is that the expected reflection continuum from 
highly ionized matter is absent in the spectral data. The 
fit with the inclusion of a reflection component from an 
ionized face-on disk shows that this component is not re- 
quired by the data, and we derive a 90% confidence upper 
limit on the subtended solid angle Jl/27r of the reflecting 
material of 0.3. Moreover, the fit is of poor quality be- 
cause reflection cannot reproduce a sharp absorption edge 
at 8.6 keV as observed (it produces smeared edges, Ross 
et al. 1999). The physical problem is that the observed 
accretion rate is too low for the required ionization level 
in the disk. In order to have He-like iron the accretion rate 
in Eddington units m — M/AfEdd must be 0.4 (> 0.2 
for >FeXXII; Matt et al. 1993). The observed bolomet- 
ric luminosity of the nucleus gives m ~ 4 x 10~* only 
(see Sect. 4.2 below). So, m fails severely the requirement 
ra > 0.2. Therefore the line and the edge are unlikely to 
be produced by reflection from a highly ionized accretion 
disk. 

We suggest instead that they are produced by trans- 
mission through a highly photoionized medium, located 
close to the nucleus. If so, this points out an interesting 
similarity with the presence of a warm absorber in Seyfert 
I's (Nandra & Pounds 1994). Usually Seyfert I's show the 
oxygen absorption edges, superimposed onto the power 



4-2. Origin of the X-ray continuum 

In the framework of standard a-disks (Shakura & Sun- 
yaev 1973) X-ray emission with a power law shape arises 
from inverse Compton scattering, by a hot optically thin 
plasma in a corona, of UV/soft X-ray photons produced by 
the accretion disk (see, e.g., Svensson 1996 for a review). 
The observed bolometric luminosity of the nucleus of M81 
corresponds to a low accretion rate. Adopting the central 
mass value of 4 x 10^ Mq, the bolometric luminosity of the 
nucleus obtained by integrating the observed spectral en- 
ergy distribution from the radio to 10 keV {L = 2.1x 
erg s"^. Ho 1999) and a radiative efficiency of 0.1, it turns 
out that m ~ 4 X lO"**. A problem with this scenario is 
that we do not detect signs of reflection from optically 
thick cold material, usually found in Seyfert I's, and at- 
tributed to the presence of an accretion disk. These signs 
are a 6.4 keV emission line with typical EW of 100-150 
eV, and a broad bump peaking at 10-20 keV (Nandra & 
Pounds 1994). As a further comparison with the contin- 
uum of classical, more luminous AGNs, note that also the 
'big blue bump', traditionally attributed to thermal ra- 
diation from an accretion disk, is absent in the spectral 
energy distribution observed for the nucleus of M81 (Ho 
1999). So, we conclude that this LLAGN is not a simple 
extension of high luminosity ones. 

An alternative possibility for the origin of the con- 
tinuum emission could be the presence of an advection 
dominated accretion flow (ADAF, Narayan & Yi 1995), 
a solution devised for AGNs radiating at very low Ed- 

^ In order to have an ionization level >FeXXII to peak, the 
ionization parameter ^>100 erg cm s~^, while ^ ~ 120 erg cm 
for He-like iron to peak. These values refer to an input 
spectrum of a power law with F = 1.5 from 13.6 eV to 13.6 
keV, passing through an optically thin shell of gas (Turner et 
al. 1992). 
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dington ratios (e.g., for NGC4258, Lasota ct al. 1996). In 
ADAF models at low mass accretion rates {m<10~^) the 
X-ray emission has a significant thermal bremsstrahlung 
contribution produced by the electrons in the flow, at 
kT ~ 10^ — 10^° K; this corresponds to a very flat power 
law (r<1.3) up to several tens of keV. As m gradually ap- 
proaches and exceeds 10"*^ the ADAF model predicts that 
inverse Compton scattering of soft synchrotron photons by 
the flow electrons will dominate the X-ray emission more 
and more. This second regime of the model is required to 
be at work for M81, in order to explain the steep spectral 
shape revealed by BeppoSAX (Sect. 3.3.1; this is possible 
when assuming a radiative efficiency < 0.1). 

In conclusion, the origin of the power law emission is 
an open question. Prom X-ray, UV and optical observa- 
tions there is no evidence for an accretion disk, and the 
conditions for an ADAF are satisfied, but the slope of the 
X-ray continuum is very similar to that of bright Seyfert 
galaxies, and it remains to be investigated whether it can 
be reproduced by an ADAF. 

4-3. Comparison with ASCA results on LINERs 

How common are the results presented here among LIN- 
ERs? On the basis of ASCA data, a canonical model was 
found to work for LLAGNs, LINERs and starburst galax- 
ies by Ptak et al. (1999). This model consists of a soft 
component of kT ^ 0.5 — 1 keV, usually spatially ex- 
tended and with absorption consistent with the Galac- 
tic value, plus a power law of mean V ~ 1.7, with small 
or no intrinsic absorption, or another thermal component 
oi kT ^ & keV. This model turned out to be success- 
ful also in the case of the detailed study of the LINERs 
NGC4579 (Terashima et al. 1998) and NGC4736 (Roberts 
et al. 1999). The results obtained here from BeppoSAX 
data are in agreement with this canonical model, with the 
addition that we can firmly establish that a power law is 
to be preferred to the hard thermal component. The pres- 
ence of a 6.7 keV gaussian line, and no evidence of a 6.4 
keV line, was also found from ASCA data for the LINERs 
NGC4579 (Terashima et al. 1998) and NGC4736 (Roberts 
et al. 1999). The EWs turned out to be larger than found 
here for M81 (360^^^^ cV, and 572^!^^? eV respectively). 

Also the absence of rapid variability on timescales less 
than a day, a likely case for the nucleus of M81 (Petre et 
al. 1993), seems to be a general property of the LINERs 
observed with ROSAT aivd ASCA, while changes in the 
X-ray luminosity over long timescales are not uncommon 
(Reichcrt et al. 1994, Ptak et al. 1998). This has been 
interpreted in terms of an ADAF being at work (Ptak et 
al. 1998). 

5. Conclusions 

The (0.1-100) keV study of M81 based on BeppoSAX 
data has revealed the presence of a power law compo- 



nent, of emission and absorption features due to highly 
ionized iron, and short term variability with a period of 
~ 2 days. These characteristics point to the presence of a 
LLAGN, as suggested by 196 based on ASCA data. Dif- 
ferences and similarities with the X-ray spectra of more 
luminous Seyfert galaxies have been investigated. 

The main result of the present work is that a 6.7 keV 
emission line, corresponding to Ka emission from He-like 
iron, is clearly detected (confirming a previous ASCA re- 
sult) together with an absorption edge at ^ 8.6 keV, due 
again to highly ionized iron. The ionization level is con- 
sistent for the line and the edge, and perhaps the same is 
true also for the column density of ionized material. So, we 
suggest that line and edge are produced by a highly ion- 
ized absorber, that could be located within 10^^ cm from 
the source of ionizing photons. Such an optically thin gas 
that transmits the X-ray photons represents a similarity 
with the warm absorber in Seyfert I's. 

Concerning the main spectral property of the contin- 
uum, BeppoSAX data reveal the presence of a single 
power law of photon index V ~ 1.84 in the (0.1-100) 
keV energy band. This is modified by cold absorption of 
A^H ~ lO'^^ cm"^, that could be entirely produced by the 
Galactic hydrogen column density plus some intrinsic A^h 
corresponding to the reddening observed for the broad line 
region. The F value, and the absence of high intrinsic ab- 
sorption, suggest that M81 might be an example of a low 
luminosity Seyfert 1. However, there are also differences 
with the X-ray properties of Seyfert I's: there is no sign of 
reflection from optically thick cold material, like the 6.4 
keV emission line and the broad bump peaking at 10-20 
keV. Even a blue bump is absent in the spectral energy 
distribution of the nucleus of M81, all this suggesting that 
this nucleus cannot be considered a simple extension to 
low luminosities of classical Seyfert I's. 

Alternatively to the standard accretion disk scenario, 
the continuum could be produced by an ADAF. This 
hypothesis is consistent with the low accretion rate in- 
ferred from the luminosity and the estimate of the central 
black hole mass; it remains to be investigated whether the 
ADAF scenario can reproduce a Seyfert-like power law 
with r ~ 1.84 over 0.1-100 keV for M81. 
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